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The time dispersion of the averaged conductance G(t) of a mesoscopic sample is calculated in the 
long time limit when t is much larger than the diffusion travelling time to- In this case the functional 
integral in the effective supersymmetric field theory is determined by the saddle point contribution. 
If t is shorter than the inverse level spacing A (At/h <C f), then G(t) decays as exp[—t/to]- In 
the ultra-long time limit (At/h 3> 1) the conductance G(t) is determined by the electron states 
that are poorly connected with the outside leads. The probability to find such a state decreases 
more slowly than any exponential funcion as t tends to infinity. It is worth mentioning, that the 
saddle point equation looks very similar to the well known Eilenberger equation in the theory of 
dirty superconductors. 

PACS numbers: 73.40. Gk, 05.45. +b, 73.20.Dx, 72.20.My 



We consider long time relaxation phenomena in a dis- 
ordered conductor that is attached to ideal leads. For 
simplicity we assume that the electrons in this conductor 
do not interact with each other, the temperature is zero 
(T = 0) and there are no inelastic processes. The total 
current I(t) at time t depends upon the voltage according 
to the Ohm law: 

I(t)= f dt'G{t-t')V(t'). 



We are interested in the asymptotic form of the conduc- 
tance G(t) as t — ► oo. 

The same problem has been considered earlier by Alt- 
shuler, Kravtsov and Lerner (AKL) Q. Our initial goal 
was to obtain their results by means of a more direct 
calculation. At present, we can neither confirm, nor dis- 
prove the AKL results. We have found an intermediate 
range of times, where the conductance G(t) decays more 
slowly than it was predicted in Ref. [El . The AKL asymp- 
tote could be valid at longer times (see discussion below). 

There are three time scales in the problem: 

1 . The mean free time r = l/vp, where vf is the Fermi 
velocity and / is the mean free path. This time scale 
determines the dispersion of the Drude conductiv- 
ity cr ~ e _t / T 

2. The time of diffusion through the sample try = 
L 2 /D, where D = 1 2 /3t is the diffusion coefficient, 
and L is the sample size. 

3. The inverse mean level spacing h/A — TivV , where 
v is the density of states and V is the volume of the 
sample. 

In a macroscopic sample the inequality r <C tu h/A is 
valid, provided that L 3> I and the disorder is weak. In- 
deed, the product Ato is connected with the dimension- 
less conductance of the sample g — 2-Kh/(tryA), which is 



large for a weak disorder. The times try and h/A enter 
into time dispersion only due to quantum corrections to 
conductivity. 

At times t <C h/A an electron can be considered a 
wave packet of many superimposed states propagating 
semi-classically. Therefore, it is natural to assume that 
the conductance G(t) is proportional to the probability of 
finding a Brownian trajectory that remains in the sample 
for the time t. For t ^ try such a probability decays as 
exp [—t/try]. Our calculations confirm this result. 

In the opposite limit, for f > h/A, the conductance 
G(t) is proportional to the probability of finding an elec- 
tron state with the life time t. In order to trap an electron 
for a long time the state must be poorly connected with 
the leads (nearly localized). We show that the probabil- 
ity of finding such a state decays non-exponcntially with 
time. Namely, G(t) ~ exp [— <?log 2 (iA)] for d = 1 and 
Git) ~ (tA)~ 9 for d — 2. These results are not valid in 
the very long time limit. We discuss this later together 
with the question of dimensional crossover. 

Instead of calculating the conductance as a function of 
time, we could have worked in the frequency representa- 
tion. In that way we would have found a singularity in 
G(ui) as uj — > 0. This singularity, however, does not affect 
the value of the d.c. conductance and therefore has an 
obscure physical meaning, while the time domain results 
have the direct interpretation. 

Since the long time asymptote corresponds to the rare 
events when the electron is nearly trapped in the sample, 
it is natural to use the saddle-point approximation. We 
carry out the following program: 

1. Express the averaged conductance ^ as a func- 
tional integral over supermatrices || (see for 
review) : 



G(t) = G e 



-t/r 



1 



/ £ e " Mt / v ^ p ^ ex P Ml ' 0) 

A=^ J dr Str{L>(Vg) 2 + 2iuuKQ} 1 



2. Vary the action A with respect to Q, taking into ac- 
count the constraint Q 2 = 1, and obtain the saddle- 
point condition which recalls the diffusion limit of 
the Eilenberger equation raj: 



2DW{QWQ) +iuj [A,Q] = 



(2) 



3. Derive the condition at the boundary with the lead 

Q I lead = A- (3) 

4. Perform the integration over u> in Eq. (Q) and ob- 
tain the self consistency condition 



y-Str(AQ) = - — 



(4) 



which allows us to exclude lu from Eq. 



5. Substitute the solution of Eq. (||) with boundary 
conditions ^ in Eq. ([!]) and obtain the results with 
exponential accuracy. 

The 8x8 supermatrix Q has commutative and anti- 
commutative matrix elements. Since Q 2 — 1 it can be 
chosen in the form |jj : 



Q = V~ l HV, V 
1 
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A = 



-1 
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This decomposition allows us to present the action A 
in the form 

A = lf J dr StT { D( y H ^ + 13 A/2 + 2iujAH }> ( 6 ) 

where M = [V _1 VV,i?]. The minimum action is 
reached for V — const, and Eq. (^) may be expressed 
in terms of 9- variables only: 



A = 



£ J dr{2 [D{V9f - 2iwcosh<9] 

-[D{X79+) 2 + 2iuucos9+} 
-[D(V6L) 2 + 2kjcos6L]} 



(7) 



where 6± = 6% ± 62- Consequently, Eq. (2) has the form: 



DV 2 + iwsinh = 0, 
DV 2 9 ± +iujsm9 ± = 



(8a) 
(8b) 



The boundary condition (||) follows from the fact that Q 
does not fluctuate in the bulk electrodes, Q = A. Hence, 
at the boundary with the ideal lead 9 = 9± = [||. 
The time decay of the conductance G(t) ~ exp(— icot) 
corresponds to real and positive values of ilu. The per- 
mitted values of frequency lu in Eq. ( pb| ) are bounded 
from below by the value lu\ ~ which corresponds 

to the linearized form of Eq. (glj). For smaller frequen- 
cies u> < lo\, which will turn out to be the only relevant 
ones, Eq. (8b) has only trivial solutions 0± = 0. Thus, 
the self-consistency equation (||) has the form: 



(9) 



The solutions of Eq. ( pa] ) depend on the sample geometry. 
We start by considering a one dimensional wire of length 
L, attached to ideal leads at x = ±L/2. If t A <C 1, 
then, to satisfy the self-consistency condition (||) we must 
choose 9 <C 1. Therefore, Eq. (|Sa| ) can be linearized. The 
solutions that satisfy the boundary conditions is 



Ccos(jmx/L), u) n 



to 



(10) 



where n is an arbitrary integer. The above formula for 
the frequency implies that in the discussed regime 



G(t) 



e-*"** = exp 



(ii) 



To obtain this result we determine the amplitude C from 
the linearized self-consistency equation, and then substi- 
tute ( |Io|) into the action A. 

For arbitrary times Eqs. ( 
coordinates have the form: 

d 2 9 7 2 . 



a) and 
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7 = 2iujtu- 



(12) 



(13) 



The solution of ( |12| ) is symmetric = 9{—z), and in 
the region 2 > is given by the quadrature: 



1 



d9' 



7 Je(z) V cosh 6*o — cosh W 



(14) 



o = 0(O) = 21ogi + 21oglogi, for 7 <l (15) 

7 7 

The function is almost linear 9 — 9${1 — 2\z\) 
everywhere except in the region \z\ < l/log(l/7) <C 1. 
Substituting Eq. Q into Eqs. (|) and (g) we get 
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T log T< 



G(t) ~ exp 



-5 log' 



(16) 



As mentioned earlier, the contribution from the indi- 
vidual nearly localized states dominates in G(t) when- 
ever tA/h 3> 1. The square modulus of the wave func- 
tion for such a state |\1/| 2 equals cosh6>. As we can 
see, this value decays exponentially towards the leads, 
where \^{x = ±L/2)\ 2 = cosh0(±l/2) = 1. Because 
of the latter condition, the current through the wire is 
equal to unity. Therefore, the escape time t is pro- 
portional to the normalization integral. This is exactly 
what is stated in the self-consistency condition (|l^) for 
1. To summarize, the wave function is localized 
in the region \x\ -C £ *C L with the localization length 
£ = L/ log(tA/h) and the probability to find such a state 
is given by Eq.(|l6|). For very long times, when £ be- 
comes less than the transverse size of the sample, the 
one-dimensional regime crosses over to a two- or three- 
dimensional one. 

In the two-dimensional case we consider a mesoscopic 
disk of radius R surrounded by a well conducting elec- 
trode. The Laplacian operator in Eq. (^]) is now two- 
dimensional and the boundary condition is 9(R) = at 
the circumference of the disk. It is natural to assume 
that the minimal action corresponds to that depends 
on the radius only and, therefore, obeys the equation: 



+ 6'/z + iu)t D smh9 = 0, 0(1) = 



(17) 



where z = r/R and to = R 2 /D. 

For t < fi/A , Eq. @ can be linearized. Its solution 
is the Bessel function 



9 = CJ (jz), 7 = V^d = Vn, 

where fi n denotes the ro-th zero of the Bessel function. 
The conductance is 

G(t) -exp > t D -^t<^h/A. (18) 

For a long time tail t 3> ?i/A, the non-linear term in 
Eq. dl7| ) is large near the origin and can be neglected 
elsewhere. As a result, 



9(z) = Clog 



1 



(19) 



for all but very small z. On the other hand, for z« 1, 
the parameter 9 is large and sinh0 = e 9 jl. With this 
approximation the solution of Eq. (|lj]) can be found j|] 
having the asymptote 

9{z) = -6>(0) + 6 log 2 + log -i + 4 log - , (20) 

7^ z 

for 7 <C z < 1. Comparing with Eq. (|l9|), we have 9(0) = 
6 log 2 + log(4/7 2 ) and C = 4. To calculate the integral 
in the self-consistency equation 



At f 1 

— -= / {cosh0(z) - \}zdz (21) 
2nh J 



we multiply Eq. (17) by z, integrate in the limits and 
1, and obtain 



dB_ 

dz 



+ iojtD / smh9(z)zdz = 0. 



(22) 



Since 9(0) 3> 1, we neglect the difference between the 
integrals in Eqs. (|2l]) and ((2^), and with asymptote jl9| ) 
finally get iu> = ig/t. The action A is dominated by the 
contribution of the tail (O) : 



A = 4. g log 



tA 

2^h' 



G 



AV 5 

At) 



(23) 



The characteristic size of the averaged 2D wave function 
is £ = jR = R(h/tA) 1 ^ 2 . The crossover to a 3D case oc- 
curs when £ becomes comparable with the film thickness. 

The consideration of the 3D case makes relevant the 
question of the validity of the diffusion approximation. 
As before, we consider a disordered drop of radius R sur- 
rounded by a well conducting lead. Analogously to what 
has been done in the 2D case, the function 9 depends on 
the radius r only and obeys Eq. (|8a|), where the Laplace 
operator is substituted by its 3D radial component. The 
boundary condition is 9(r = R) = 0. The analysis of the 
linear regime is similar to that for ID and 2D cases and 
gives for t D = R 2 /D < t < fi/A: 



A = ir 2 t/t D , G(t) ~ exp(— 7r t/to) 
The nonlinear in 9 regime leads to the equation 



dz 2 



2 d9 

z dz 



iutu sinh0 = 0, 



r 
R 



(24) 



(25) 



The analysis of this equation shows that the permitted 
values of uj are larger than a certain value coq > 0, and 
that the integral in Eq. (^ remains finite even for the 
solutions of Eq. ( p5| ) with 9(r = 0) — > oo. As a result, 
the self-consistency equation cannot be satisfied for suf- 
ficiently long time t > ft/ A. Thus, for \ojto\ <§C 1, all 
non-trivial solutions of Eq. ( |25| ) satisfying the condition 
0(1) = are singular at z — * 0. Therefore, the derivative 
d9/dr becomes comparable with the inverse mean- free 
path 1 jl for a certain radius r* . The diffusion approxima- 
tion inevitably breaks down for smaller distances, where 
non-local corrections become important. 

We do not try to solve the kinetic problem now but as- 
sume that the mentioned non-locality smoothes out the 
singularity at the origin. We also assume that, similarly 
to what has happened in the ID and 2D cases, the non- 
linear term in Eq. ( p5| ) can be neglected at r > r* and is 
important for r — r*. Thus 
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0(r) ~ C(- - 1), r > r, 
r 

d6» _ OR 
- ^ | r=P . - -jj-, 

and r, = {CIR) 1 ' 2 . Then 0, = 0(r„) = {CR/l) 1 ' 2 and, 
finally, 



wt_D exp(0») 



t R 2 1 /i? 



which gives 
= lo, 



1 /i? 



1 fR 



Using the self-consistency condition we express the fre- 
quency lj through the time t and obtain a rough estimate 
for the action 



log Q 



G(t) ~ exp [-A{t)\ , (26) 



where A = h/pp is the Fermi wavelength and the expo- 
nent a ~ 1 can only be determined from the solution of 
the kinetic problem. Equation (Eq) looks similar to the 
AKL result if a = 2. The action (Pqj becomes compara- 
ble with that from Eq. (|3) at 



t — — tjj 



h l 
AR' 



t D < U < 



A 



(27) 




Therefore, we have two independent contributions to 
the conductance G(t), which is dominated by one of them 
( Eq. @) at t < U and by the other (Eq. @) at * > i*. 

It is sensible now to analyze whether the diffusion 
treatment is valid in the ID and 2D cases. Using the 
solutions of Eqs. (|l^) and (|l7|) we find the value of 
such that for t < the derivative ld6/dr is less than 
unity. This gives: 



(28) 



Therefore, we can expect that the asymptotes (|T(|) and 
( p3| ) are valid for t < t*. At longer times t > i» for all 
dimensions d = 1,2,3 the asymptote cannot be found 
within the diffusion approximation. A detailed kinetic 
analysis of this problem will be done in a separate work. 

We want now to speculate on the possible relation of 
this long time tail to the AKL asymptote. AKL studied 
the coefficient growth rate in a power expansion of G(uj) 
in lot. The rate they found means that G(t) has a univer- 
sal logarithmically-normal asymptote at very long times. 
The presence of the mean free time r in this expansion 
makes it plausible that the AKL asymptote is related to 
a kind of kinetic problem. This is why we hope that at 
t> the tail will match the results of Ref . |l[] . 



In the case of a tunnel barrier at the sample-lead in- 
terface , the time dispersion of the conductance can be 
considered in the same way with the usage of the gener- 
alized boundary condition (see Ref. H). 
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The substitutions ; 
Eq. ([H]) into the form tp 
integral. 



e 17 and 6 — ^(r/) + 2r/ transform 



+ 7 2 /4e* = which has the first 
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